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Abstract 

Imposing non-integrable constraints on Ricci flows of (pseudo) Rie- 
mannian metrics wc model mutual transforms to, and from, non- 
Riemannian spaces. Such evolutions of geometries and physical the- 
ories can be modelled for nonholonomic manifolds and vector/ tan- 
gent bundles enabled with fundamental geometric objects determin- 
ing Lagrange-Finsler and/or Einstein spaces. Prescribing correspond- 
ing classes of generating functions, we construct different types of 
stochastic, fractional, nonholonomic etc models of evolution for non- 
linear dynamical systems, exact solutions of Einstein equations and/or 
Lagrange-Finsler configurations. The main result of this paper consists 
in a proof of existence of unique and positive solutions of nonlinear dif- 
fusion equations which can be related to stochastic solutions in gravity 
and Ricci flow theory. This allows us to formulate stochastic modi- 
fications of Perelman's functionals and prove the main theorems for 
stochastic Ricci flow evolution. We show that nonholonomic Ricci flow 
diffusion can be with self-organized critical behavior, for gravitational 
and Lagrange-Finsler systems, and that a statistical/ thermodynamic 
analogy to stochastic geometric evolution can be formulated. 
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1 Introduction 

The purpose of this paper is to analyze models of Ricci flows of random 
metrics with self-organize criticality in general relativity (and various mod- 
ifications, for instance, of Lagrange-Finsler type theories of gravity) in the 
framework of stochastic evolution equations. In other words, we shall ana- 
lyze scenarios of stochastic/diffusion Ricci evolution of curved spaces with 
rich geometric/physical structural resulting in random gravitational systems 
with critical points as attractors. 

This is a partner work to [T]; our approach to stochastic solutions and dif- 
fusion of gravitational fields is elaborated in the framework of the theory of 
stochastic evolution equations following the methods elaborated in j2j[3j|4l[5] 
(see also references therein) but modified/adapted for curved spaces and 
nonholonomic Ricci evolution of Riemannian and non-Riemannian geome- 
tries, PHI EE]). 

1 when, for instance, metrics are solutions of Einstein equations with certain prescribed 
symmetries and/or nonholonomic constraints, solitonic configurations, stochastic proper- 
ties etc 
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The theory of stochastic processes and diffusion on curved spaces has 
been studied in mathematics and physics from different perspectives re- 
lated to new directions in differential geometry and partial differential equa- 
tions, geometric analysis and evolution theory, kinetic and thermodynamic 
processes with local anisotropy, and various applications in cosmology and 
astrophysics. We refer to a series of works containing original key ideas, 
methods and reviews of results P33 [III Q21 1131 D31 H51 HSl HZl HBl HSl EDI EU - 

In our recent paper [I], we studied exact solutions in gravity defined by 
stochastic generating functions and nonholonomic diffusion of gravitational 
fields. Such generic off-diagonal metrics (which can not be diagonalized 
by coordinate transforms) describe generic nonlinear and random processes 
with very complex structure of gravitational "spacetime ether", for instance, 
solitonic interactions \22\ [23] , fractional configurations |24[ [25| [26] and non- 
linear wave interactions and geometric evolution |27|, [2"51 [29] . 

In this work, the gravitational stochastic Ricci flow and field interactions 
phenomena with self-organized criticality are studied following a synthesis of 
the method of anholonomic deformation/frames of constructing exact solu- 
tions in gravity [301 ED E21 [331 EI] an d the theory of nonholonomic diffusion 
on curved spaces (additionally to above mentioned works, we cite our con- 
tributions on nonholonomic diffusion in (super) vector bundles and locally 
anisotropic kinetics and thermodynamics [131 E3 O EH] ) . To the best our 
knowledge this is the first attempt to provide an unified geometric formalism 
for locally anisotropic diffusion processes and Ricci flows of Einstein/-Finsler 
spaces, when certain classes of nonholonomic constraints result in stochastic 
evolution equations which can be approached in a mathematically strict way 

[21E1S1E]! 

This paper is organized as follow. In section [2] we provide an introduction 
to the theory of Ricci flow evolution of nonholonomic geometries. Section 

2 During last two decades, self-organized criticality is widely studied in physics 
[361 1371 1381 1401 I41j . Our idea is that such effects are possible for gravitational fields 
and their evolution being derived as corresponding solutions. Classical vacuum structure 
in general relativity may have various sophisticate nontrivial topological and geomet- 
ric configurations with possible stochastic, solitonic, instanton, black hole, wormhole etc 
properties. We can say that a curved spacetime is modeled by a nontrivial fundamental 
gravitational "ether" which, as a matter of principe, may have a very complex nonlinear 
sure and/or random behavior which may result in configurations with self-organized crit- 
icality. The term "porous" media is a conventional one for gravitational configurations 
related to corresponding equations which are formally similar to certain equations for real 
porous matter. In this work, we show in explicit form what type of nonholonomic con- 
straints we have to impose on gravitational field equations, and their possible geometric 
evolution, in order to have critical points as attractors. 
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[3] is a development of some recent results on stochastic diffusion equations 
(in our case) on nonholonomic manifolds and to proof of existence of unique 
and positive solutions for such systems in gravity and geometric mechanics. 
In section 0] we develop the theory of stochastic nonholonomic Ricci flows: 
there are considered stochastic modifications of Perelman's functionals and 
proven the main theorems for stochastic evolution equations. Finally, a 
statistical analogy for stochastic Ricci flows is proposed. 

Acknowledgement: The author is grateful to Academician V. Barbu 
for the opportunity to present a part of results of this paper at the Seminar of 
Institute of Mathematics "O. Mayer' of Romanian Academy, Ia§i, Romania. 

2 Ricci Flow Evolution and Nonholonomic Mani- 
folds 

The Ricci flow theory [421 l4"3j [4"4"1 [i"5] , related geometric analysis and var- 
ious applications (see [Ml 071 05] for reviews of results and methods) became 
one of the most intensively developing branch of modern mathematics. The 
most important achievement of this theory was the proof of W. Thurston's 
Geometrization Conjecture by Grisha Perelman |43 t 144 1 ^5]- The main re- 
sults on Ricci flow evolution were proved originally for (pseudo) Riemannian 
and Kahler geometries. In a series of our works [6l[3[8j|9l[2HI39j[50], we stud- 
ied Ricci flow evolutions of geometries and physical models (of gravity with 
symmetric and nonsymmetric metrics, geometric mechanics, fractional and 
noncommutative generalizations) when the field equations are subjected to 
nonholonomic constraints and the field/evolution solutions, mutually trans- 
form as Riemann and generalized non-Riemann geometries. In this section 
we outline the anholonomic deformation method of constructing exact solu- 
tions in gravity and Ricci flow theory in a form necessary for further devel- 
opments (in next sections) for stochastic evolution and Ricci flow diffusion. 

2.1 Modelling Einstein and Finsler geometries on nonholo- 
nomic manifolds 

In the following we shall consider a nonholonomic manifold V = (V,Af) 
enabled with nonlinear connection (N-connection) structure N : TV — > 
hV © t>V as a nonholonomic distribution J\f = N defining a conventional 
horizontal (h) and vertical (v) splitting on V. We consider that V is a four 
dimensional, 4-d, (pseudo) Riemann manifold of signature (+,+,—+). For 
models of Lagrange-Finsler geometry, we can take V = TM to be the total 
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space of a manifold M, when the typical fiber in such a manifold is provided 
with a hyperbolic structure in order to mimic a local (pseudo) Euclidian 
structure. We shall use the notations from [I] (readers may find details on 
the geometry of nonholonomic manifolds [511 [52"1 I53| ) in monograph |54| ; 
for purposes of modern geometry, mechanics and physics and mathematical 
relativity, the formalism is developed in [561 E21 EQ]; the concept of N- 
connection is contained in coefficient form in [55] being developed in details 
on a number of works on Lagrange-Finsler geometry and generalizations/ 
modifications [571 EH E3 El E]). 

We label local coordinates on V in the form u a = (x l ,y a ), were i, j, ... = 
1, 2 and a, 6, ... = 3, 4 (in brief, we write u = (x, y)); similar values are taken 
by any variants of primed, underlined etc indices, for instance, a = (i,a) 
and j3' = (j',b'). A local frame and co-frame are written e a = (ej,e a ) and 

= (e', e b ), when frame transforms are parametrized e a = e a a (u)e a >. 

2.1.1 Nonlinear connections and Einstein equations 

A N-connection structure on V, with local coefficients {N®} stated with 
respect to a coordinate basis, allows us to define the so-called N-elongated 
(equivalently, N-adapted) frames, i.e. , respectively, partial derivatives and 
differentials, 

= d i -N?d a ,e b = d b = -^\ , (1) 

= [j = dx\e a = dy a + Nfdx 1 ) . (2) 

With respect to such bases, the geometric objects are written in N-adapted 
form and called as distinguished objects (in brief, d-objects), for instance, 
d-vectors, d-tensors, d-connections etc. 

Our geometric arena consists from nonholonomic manifolds / bundles (we 
shall use as equivalent the terms, spaces, or spacetimes, for corresponding 
signatures) given by data (N, g, D) , where the d-metric g is parametrized 
in the form 

g = g ij dx i dx j + h ab (dy a + N%dx k )®(dy b + N%dx k ), (3) 



and D may be taken to be the canonical d-connection D = yhD, vD 

uniquely defined from the conditions that it is metric compatible, Dg = 0, 
and with vanishing "pure" h- and v-components of torsion T of D. Even 
there are nontrivial h-v components of torsion, there is a distortion relation 

D = V + Z, (4) 
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where all values, i.e. the canonical d-connection, D = {r°jg }, the Levi- 
Civita connection V = {r®^ } (subjected to the conditions that it is torsion- 
less and Vg = 0) and the distortion, Z = {Z°g ~ T"^^}, are completely 
and uniquely defined by the same metric tensor g. We can work equivalently 
with both linear connections D and/or V (the first one is N-adapted but 
the second one is not)H 

The Einstein equations in general relativity can be written for D, 

1 

R- ps ~ 2§/3<5 S ^ = T /3<5> ( 5 ) 
L c aj = e a (N*), q 6 = 0, f^, = 0, (6) 

where R ps is the Ricci tensor for T' y aj3 , S R = g^R p§ and the source d- 

tensor ¥ps is constructed for the same metric but with D (formulas are 
similar to those in general relativity with V; we consider ¥ p$ — > nTp$ for 
D — > V). If the constraints ([6]) are satisfied the distortion d-tensors Z 1 ^ 

from (jl]), determined by d-torsion T 7 Q ^, became zero and ([5]) are equivalent 
to the "standard" equations 

R I3S- T^gpsR = xTps, (7) 

written for the Levi-Civita connection V = {r 7 Q ^}. In formulas ([7]), R ps 
and R are respectively the Ricci tensor and scalar curvature of V; it is also 
considered the energy-momentum tensor for matter, T a p, where x = const. 

Here we note that it is not possible to integrate analytically, in general 
form, the Levi-Civita form of the system of partial differential equations (|7|) 
because of it generic nonlinearity and complexity. Nevertheless, it is possi- 
ble to solve in very general forms the version for the canonical d-connection 
(and for the Cartan d-connection in Finsler gravity) the system ([5]), see de- 
tails in Refs. [331 EH EI], which motivates the idea to introduce the so-called 
Lagrange-Finsler variables in general relativity. Imposing additionally the 
conditions ©, we can restrict the nonholonomic integral varieties and gen- 
erate solution of Einstein equations for the Levi-Civita connection. 

2.1.2 Lagrange Finsler structures and Einstein gravity 

Let us consider how data (N, g, D) can be parametrized on V = TM 
in such a form that they will define a Lagrange and/or Finsler geometry. 

3 The formulas for N-adapted coefficients of geometric objects (|4j are given in many our 
works (see, for instance, [Tl 1321 [30] and [56], for explicit constructions in Lagrange-Finsler 
gravity). For simplicity, we omit such considerations in this work. 
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We introduce frame transforms L g a '/3'('u) = e a a ,(u)e «(u)g a/ g(u), when the 
coefficients of d-metric ([3j) are transformed into coefficients of 

L g= L g ij dx i ®dxP + L h ab (dy a + L N%dx k )®(dy b + L N b k dx k ), (8) 

where L gij ~ L /i a 6 = \ & QyaQ y \ ^ an d L iV^ are respectively the Hessian and 
canonical N-connection defined by a regular Lagrangian L(u) = L(x r ,y c ), 
see details in [561 [32], [30] (and, for Ricci flows of Lagrange-Finsler struc- 
tures, [8]). The d-metric ([8]) is called the Sasaki lift on TM of a regular 
Lagrange structure. In arbitrary frames/coordinates, a Lagrange d-metric 
L g is represented in the form g ([3]) and, inversely, any g can be transformed 
into a ^g, into certain Lagrange variables for a correspondingly chosen dis- 
tribution/generating Lagrange function L(u). The formal constructions can 
be performed geometrically in similar forms on (pseudo) Riemannian space- 
times and tangent bundles but have different physical interpretations (in the 
first case, we can elaborate analogous Lagrange-Finsler models of Einstein 
gravity but in the second case it is provided a Finsler similar geometrization 
of regular Lagrange mechanics). 

We note that taking L = F 2 (x, y), where the homogeneous on y-variables 
real function F is the fundamental/generating Finsler function (see rigor- 
ous mathematical definitions and details in [561 1571 158] ). we model a Finsler 
geometry on a (pseudo) Riemannian manifold or on V =TM. Using the ge- 
ometric formalism of nonholonomic distributions, associated N-connections 
and the geometry of nonholonomic manifolds, we work in a unified form with 
all types of (pseudo) Riemann and Lagrange-Finsler geometries. We distin- 
guish such constructions by additional suppositions on the structure of man- 
ifold/bundle spaces and the type of linear connection, D, V or (the Cartan- 
Finsler d-connection) C D, we chose for the geometric data (N,g, D) . 

Finally, we shall say that there are used nonholonomic (Lagrange, or 
Finsler, variables) in Einstein gravity on V if the theory is defined by 

any data ~ L N; g ~ L g; D ~ D, or C D^ , which can be equivalently 

transformed, via distortion relations of type (j4|), into data (g, V) . For theo- 
ries on TM, we can elaborate Einstein-Lagrange /-Finsler like theories (with 
metrics and connections depending on velocity type coordinates, y a ) using 
instead of V any convenient for physical purposes D, or C D (or any other 
metric compatible d-connection completely defined by g and adapted to 
a chosen N). In abstract form, a N-adapted Einstein model on a general 
nonholonomic manifold V generates for corresponding nonholonomic distri- 
butions an Einstein and/or Lagrange-Finsler, spacetime geometry. 
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2.2 Nonholonomic Ricci flows of Einstein/— Finsler spaces 



The original Ricci flow theory was proposed |42j with an evolution equa- 
tion for a set of Riemannian metrics g a /3(x) an d corresponding Ricci tensors 
R a p{x) parametrized by a real parameter x- We can write the Hamilton's 
equations in the so-called normalized form, see details in [4"6 l |4"? 1 |4"8], 

2r 

Q-9aJ3 = -2 RgJ5_ + —Rgafh (9) 

describing (holonomic) Ricci flows with respect to a coordinate base d a = 
d/du—; the normalizing factor r = J RdVol/Vol is introduced in order to 
preserve the volume Vol; Ra/3 and R = g-^- Rgfi are computed for the 
Levi-Civita connection V0 Grisha Perelman's fundamental results |43l |4"4"1 
[45] were based on original idea to prove that the Ricci flow is not only 
a gradient flow but also can be defined as a dynamical system on the 
spaces of Riemannian metrics. He introduced two Lyapunov type functionals 
and proved that evolution equations of type ([9]) can be derived following a 
corresponding variational calculus. 

In our approach, we studied modified Ricci flow evolution equations, 
and the corresponding N-adapted functionals, for generalized commuta- 
tive and noncommutative geometries (Riemann and Lagrange-Finsler ones) 
[HI El El [23] . We also were interested to study Ricci flows of exact so- 
lutions in various theories of gravity [491 150] 128"! 129] . In this work, we 
shall elaborate a theory of Ricci evolution with self-organized criticality for 
stochastic solutions for Einstein spaces following geometric using the canon- 
ical d-connection D. 

If V — > D, we have to change R a p — > R Q( g which transforms equations 
([9]) into an N-adapted system of evolution equations for Ricci flows of (in 
this work, we consider symmetric metrics, see details in [6~1 [7]), 











— 9ij = 2 LA/fiVj (£ a6 - \h a b ) - R tJ + A^j - ha— (NfNf),(10) 
Kb = -2 (R ab -Xh ab ), (11) 



dx 

R ia = Oand R ai = 0, (12) 

where the Ricci coefficients R^ and R ab are computed with respect to coor- 
dinate coframes and the cosmological constant A includes the normalization 



4 We underline the indices with respect to the coordinate bases but not with respect to 
some 'N-elongated' local bases. 
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factor Tp On nonholonomic manifolds, we can prescribe any convenient for 
our purposes nonholonomic distributions. For simplicity, we shall fix such 
distributions when -j^(Nf) = and model stochastic and/or Ricci flow evo- 
lution for nonholonomic Einstein manifolds with R a f, = Xh a f, and Rij = Xgij 
for so-called "stationary configurations", when the equations © with the 
source determined by A are satisfied. The h- and v-components of Ricci 
flow evolutions of gravitational configurations for J^(N£) = , are similar to 
those derived from Perelman's nonholonomic functionals after corresponding 
nonholonomic deformationsll Such equation can be written in the form 



dgjj _ g dhab 
OX ox 



^■^iji a " — — 2i? a t,, (13) 



R ia = and R ai = 0. 

For simplicity, we shall analyze Ricci flows of families of ansatz for d- 
metrics parametrized in the form 

x g = e^dx* ®dx { + h 3 {x k ,t, X )e 3 ®e 3 + h 4 (x k ,t, X )e A ®e' i , (14) 
e 3 = dt + w i (x k ,t,x)dx t ,e' i = dy 4 + rii(x k ,t,x)dx\ 

The system of equations (fT3j) for R^ = and T2(x k ,v) = A and existing 
one Killing symmetry, on vector d/dy 4 (the coefficients do not depend on 
variable y 4 ), with evolution only of the v-parts, transform into 

^ + V" = 0, (15) 

a~ h s = — ,^-/i4 = 7—, (16) 

OX h 4 Ox h 3 

(3 Wi + oii = 0, (17) 

nf + 7< = (18) 



where 



0( X ) = 4>(x k ,t, X )=ln\^=\, (19) 

V I "3«4 1 



a; 



hld i( j>, p = hX<f>*,y= (ln|/i 4 | 3/2 /N 



In the above formulas we wrote the partial derivatives in the form a' = 
da/dx 1 , a' = da/dx 2 , a* = da/dt and we shall use d x a = da/dx- The 
conditions of zero torsion, i. e. constraints ([6]), are written in the form 

w* = ej In \h 4 \,e k Wi = eiW k , n* = 0, din k = d k rn. (20) 



5 In next section, we derive such equations for stochastic Ricci flows and related Perel- 
man's functionals. 
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We have to impose additionally (|2Q|) if we wont to consider nonholonomic 
Ricci flows determined by equations f)15|) — (|X8[) in a form when the configu- 
rations for the Levi-Civita connections V(x) are extracted. 

If h\(x) ^ 0; T2 = A ^ 0, we get </>*(x) an d we can generate 
families of exact solutions of (0), with diagonal nontrivial source = 
diag[\, A, 0, 0], if 

^(x) = 2MxMx)A/0*(x). (21) 

We conclude that the family of ansatz for d-metrics (fT4"|) subjected to 
the conditions (|15p - (|21|) define Ricci flow evolutions on a real parameter 
X of a class of generic off-diagonal solutions determining correspondingly 
(non) holonomic Einstein manifolds. They depend explicitly on the type 
of families of generating functions </>(x)- If such functions are random ones 
subjected to the conditions to solve certain stochastic/diffusion equations, 
we can say that we generated stochastic Einstein equations (for details, see 
our partner work [TJ where there are also analyzed classes of solutions of 
Einstein equations with /i| = 0, or /ig = 0; the length of this paper does not 
allow us to analyze the Ricci flow evolution of such more special classes of 
Einstein manifolds) evolving, in general, randomly, on parameter x- 

3 Nonholonomic Diffusion and Self— Organized Cri- 
ticality 

We studied various examples with physically important solutions (black 
holes/ellipsoid, wormholes, pp- and/or solitonic waves etc) evolving under 
nonholonomic Ricci flows \49\ \50\ [25| [25] . Those solutions were with sure 
coefficients and sources. To our knowledge, it was not yet analyzed the 
evolution of, in general, non-Riemann geometries when certain coefficients 
of metrics and connections are random ones, i.e. the Ricci flow theory is 
with stochastic evolution. 

In this work, we shall develop such a theory of stochastic Ricci flows, for 
simplicity, for families of metrics (|14p . which allows us to generate solutions 
in explicit form and to prove the phenomena of self-organizing criticality of 
gravitational fields for various stationary and stochastic Ricci flow evolution. 
The results from papers [2j 01 [5] are crucial in proving the existence 
of unique and positive solutions. It should be emphasized here that the 
formalism of nonholonomic distributions with associated N-connections is a 
very important geometric tool for connecting the theory of Ricci flows to 
nonlinear diffusion and gravity theories and possible applications. 
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3.1 Stochastic diffusion equations on nonholonomic mani- 
folds 

We studied some examples of nonlinear stochastic diffusion equations on 
nonholonomic V in the partner work [IJ, for the so-called N-adapted (L, A)- 
diffusion and stochastic solutions of Einstein equations. It was constructed 
the corresponding Laplace-Beltrami operator for the canonical d-connection 
D. We can generate stochastic d-metrics of type (|14|) if we take random 
generating functions 4>(x) (|19p . with stochastic evolution parameter \- Using 
(|2ip . various families of exact solutions of Einstein equations can be defined. 
An important mathematical/physical problem is to state certain general 
conditions when the Ricci flow evolution equation have unique and positive 
solutions. 

Let us introduce of nonholonomic geometric and N-adapted stochastic 
calculus framework. We consider an open bounder domain U C vjf| with the 
spatial dimV < 3 with smooth boundary dU, when the Laplace-Beltrami 
operator A is completely defined by data ^N, g, . Out stochastic non- 
holonomic Ricci flows will be modelled using a nonlinear evolution equation 

5U( X ) - A*(U(x))5x 9 (r(U(x))6W(x), on (0, oo) x U, 

#(£%)) 3 0, on (0, oo) x aW, (22) 
U (0, u) = u on U. 

In the above formulas, 5W(x) is a Wiener process, the initial datum u is 
given for "rolling" the stochastic process on nonholonomic curved manifold 
V, locally on carts of a covering atlas and we can consider any maximal 
monotone (possible multivalued) graph with polynomial growth of "coer- 
cive" function ^ : M — >• 2 K . It is also possible to introduce a correspondingly 
parametrized random forcing term 

oo 

a(U)dW = J2 v kU(l,e k ) 2 e k (23) 

k=l 

for any I G L 2 (U), where (•, -fc)2 and e k are respectively the scalar product 
and an orthonormal basis in L 2 (U) (which is determined, in our case, for a 
given d-metric structure on V), when v k is a sequence of positive number 
and the set f3 k = (h e k) can be associated to a sequence of independent stan- 
dard Brownian motions on a filtered probability space ($7, J 7 , {J r x } x >o, P) . 

6 we suppose that our nonholonomic manifold is covered by such open regions 
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The physical meaning of equations (|22p depends on the type of models 
we are going to elaborate. For instance, such an equation described the 
dynamics of flows in porous media; for more general assumptions, it models 
phase transitions with melting and solidification processes in the presence of 
a random force, i. e. term. Our proposal [T] was to introduce random gener- 
ating functions for solutions in gravity (with the possibility to include noise 
for matter energy-momentum ) considering vacuum non-vacuum spacetime 
configurations as a complex "ether" media (following Mach ideas on inertial 
forces and gravity) with possible nonholonomic structure, singularities and 
diffusion of gravitational and matter fields. 

In this work, we show that gravitational fields (curved spacetime) may 
self-organize itself due to stochastic gravitational effects and /or under non- 
holonomic stochastic Ricci flow evolution. It is possible to formulate stochas- 
tic analogs of gravitational Ricci flow evolution equations (a geometric non- 
linear diffusion driven by with stochastic components of the Ricci tensor, 
with certain limits to nonholonomic versions Laplace-Beltrami operators). 

Our goal is to prove (using a synthesis of methods and results from [21 
El IU [5] and Perelman's functional approach [331 SU 05] generalized for non- 
holonomic geometries [61 [7]) that a class of stochastic evolution/gravitational 
field equations can uniquely solved in a form satisfying positive conditions. 
There will be used sure evolution analogs of equations (|22p for a function 
f(x , x 2 , v, x) which is modelled by 



OX 



DC/ 



R-S, (24) 



where R + S is the scalar curvature of D. Considering random generating 

functions we induce a gravitational random forcing term in formula (j23|) . 

2 

We say that ^ is such way defined that the terms —AU+ DC/ are analogs 

of A\£ in (|22|) when a unique solution of this equation is related to a unique 
solution of (|24p in direct form or computing certain expectation values. We 
shall also consider unique stochastic solution for any well defined functions 
f(U) = f(u) = f(x\x 2 ,v, X ). 

3.2 Existence of unique and positive solutions; physical set- 
ting for gravity and Lagrange— Finsler spaces 

On flat three dimensional real spaces, the existence problem for stochas- 
tic equations ([22]) with additive and multiplicative noise was studied in Ref. 
[3J and, with generalized conditions, in [21 0]. We shall search for general 
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conditions of existence of equations ([22]) and (p^|) on U C V under such 
assumptions: 

Hypothesis 3.1 1. The partition with IA C V and localization of oper- 
ators A (the domain of this operator is H 2 (IA) n Hq(U)), see below 
Notation \3.1\ and D are such way parametrized via nonholonomic dis- 
tributions that the coercive function is a maximal monotone multi- 
valent function from M into R, when G ^(0); 

2. 3 C > and a > u>/ien V£i w can wi£e sup{|#[ : 8 € ^(r)} < 
C(l + \r\ a ); 

3. it is possible to fix the nonholonomic distributions, the canonical d- 

^ oo 

connection D and sequence in such a form that locally ^ v k^k < - 

k=l 

+oo, for Afc being the eigenvalues of the Laplace-Beltrami d-operator 
— A on IA C V with Dirichle boundary conditions on dlA. 

Roughly speaking, the conditions of above hypothesis are defined to 
"roll" on IA, in N-adapted form, for the canonical d-connection D and A, the 
assumptions from Hypothesis 1.1 in Ref. [2]. Similarly, the equations (j22j) 
and further developments in this sections are "nonholonomic modifications" 
of the results on nonlinear diffusion from the mentioned papers [21 EJ H] but 
related to stochastic Ricci flows via profs of existence for (|24p PI 

Notation 3.1 1. We have o-{u) £ L 2 (L 2 (IA), H~ l (U)) with correspond- 
ing spaces/ conditions for all Hilbert-Schmidt d-operators from L 2 (IA) 
into H~ l (JA), when there is the Lipschitz continuity from H~ l (U) into 
L 2 (L 2 (lA),H~\lA)). 

2. By L p (U),p > 1, we denote the space of p-integrable functions with 
norm \ ■ \ p . The spaces H k ilA) C L 2 (IA), (k = 1,2), are the standard 
Sobolev spaces onU. We consider that H^ilA) is the subspace of H 1 ^) 
with zero trace on the boundary. 

3. Denoting rl as a Hilbert space, for p,q 6 [1, +00], we write 

Lw [(0, T x) ! LP f or the space of all q-integrable processes z : 

7 We recommend readers to consult those papers on nonlinear diffusion and self- 
organized criticality where basic concepts and veri important results are stated following 
a rigorous mathematical formalism. The limits of this work does not allow us to re- 
peat and develop in a detailed form those results for h- and v-splitting, with N-adapted 
constructions, on U C V. 
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[0, T x] L p (^l;T-L) are adapted to the filtration {F x } x >o- It is also 
used Cw [(O, T x) j L 2 (Q] H)] for the space of all %-valued N-adapted 
processes being mean square continuous. L(H) denotes the space of 
bounded linear operators equipped with the above introduced norm. 

4- We model 1-1 as a distribution space H = = (Hl(U)) provided 

with a N-adapted scalar product and norm defined by d-metric g and 
canonical Laplace-Beltrami d-operator, A = —A, 

( V, 2 f) = [ a- 1 7(0 2 f(0V\mm for |/u = vTJTT). 

JU 

Under above assumptions, on any hi C V and for a parameter Xi we can 
use the main result from [2] that if it € L p (U),p > max{2a,4}, then there is 
a unique strong solution to equation (|22l) (any such solution is nonnegative 
if the initial data u are also nonnegative). 

In our nonholonomic setting for gravity and flows of geometries, we can 
solve different physical problems: We do not model usual porous media, but 
a gravitational " ether" interacting as a solution of Einstein equations and / or 
following a scenarios of Ricci flow stochastic evolution. For instance, there 
are some important examples of modeling for gravitational configurations 
by nonholonomic diffusion equations of type fj22[) : 

Example 3.1 1. For certain constants °x, p, 1 a, 2 a £ (0, +oo), we 
can model a gravitational thermo-field theory with heat conduction, or 
phase transitions in spacetime "porous foam" if 

{ 1 ol(x - °x)> f or X < °x; 
[o, P ], for x = °x; 
2 «(X - °X) + P> f° r X > °X- 

2. In a spacetime with black holes and singularities, it is important to 
consider the nonlinear singular nonholonomic diffusion equation 

5U( X ) - pdiv [ °6 (£%)) DC%)] d X = a (U( X )) 5W(x), 

for N-adapted divergence div, where °5 is the Dirac measure concen- 



trated at the origin, which is possible for Vl/(x) 



PX/\x\, if X + 0; 
[-1,1], if X = 0. 



3. Choosing ^(x) = |x| as ^<7 n X with < a < 1, for a local diffusion 
problem with free boundary and a random forcing term proportional to 
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U (x) — c Xi when the value c \ determine the critical point of diffusion 
process U{x), we get the a particular nonholonomic diffusion equation 



5U(x) - A (Hev( X ) + x) [U( X ) - c u]d X = a (U( X ) - c u ) 5W( X )- 

In this formula, it is used the Heavside step function Hev(x) = {0, if 
X < 0; [0, 1], if x = 0; 1) if X > 0}- Certain physical models for self- 
organized criticality and their rigorous mathematical study are given 
in Refs. \31\ \3b\ For spacetime evolution, such criticality can be 
related with behavior of certain effective physical constants and/or dif- 
fusion gravitational phase transitions subjected to nonholonomic con- 
straints. One might consider scenarios when the supercritical region 
{U(x) > Cu } i- s absorbed asymptotically during evolution by the criti- 
cal one {U(x) = c u} . Such models seem to have applications in mod- 
ern quantum gravity (as alternative to Horava-Lifshitz phase transi- 
tions) and or in cosmology driven by Ricci flow models. 

For the above given examples, we can not apply the general existence 
theory of infinite dimensional stochastic equations in Hilbert space (enabled 
with nonlinear maximal monotone operators). We have to elaborate a INT- 
adapted approach. 



Definition 3.1 We call a solution of \22\) any H-valued continuous J~y~ 
adapted and N-adapted process U(x) = U(x,z), for z G % , on [0, T x] 
if U G L p (nx (0, T x) x Uj) n I? (0, T X ; L 2 (fl; H)) ,p > a, and 3r, G 
U>l a (n x (0, T x) x U)) such that F-a.s. 

(U(x,z), ej) 2 = (z, ej} 2 + 

j I n^OAejiOVWOl^ds + J2"k J (U(s,z) e k ,e,) 2 d^ k (s), 



u 



k=l 



Vj G N and rj G a.e. in ft x (0, 1 x) x U. 

Finally, we formulate a N-adapted generalization of the existence theo- 
rem (Main Result of [2]): 

Theorem 3.1 For each z G L p (U),p > max{2a,4} and conditions of Hy- 
pothesis^ there is a unique solution U [(0, T x) ;LP(Q;U)] , see No- 
tation HOI to \22\) . If additionally z is nonnegative a.e. in U then P-o.s. 
U( X ,z)(0>0,fora.e. (x, £) G (0, oo) x U. 
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Proof is similar to that provided for flat spaces. We have to "roll" on 
atlas carts and dub the constructions for h— and ^-components and using the 
Laplace-Beltrami operator for A determined by the canonical d-connection 
D. We omit such technical results. The most important consequence of 
this Theorem is that the existence of a unique solution for such diffusion 
equations can be related to stochastic Ricci flows. This allows us to derive in 
a unique form the evolution equations and related (stochastic) fundamental 
functionals. 



4 Stochastic Nonholonomic Ricci Flows 

For a general random generating function cp(x) introduced into a d- 
metric (|14p . it is not clear how to define the Perelman functionals and 
derive the Hamilton evolution equations for stochastic Ricci flows. The 
goal of this section is to sketch in brief a self-consistent stochastic version 
of Hamilton-Perelman theory for nonlinear diffusion of metric coefficients 
when the assumptions of Hypothesis [1] and conditions of Theorem 13.11 are 
satisfied for certain generating/normalizing functions. 

4.1 A stochastic modification of Perelman's functionals 

The Perelman's functionals were introduced for Ricci flows of Rieman- 
nian metrics and Levi-Civita connection 143, 44, 45] are written in the form 



,^(g,/) = J (<R+\Vff)e-f dV, (25) 
v 

,W(g, /, r) = j [t ( ,i? + | V/|) 2 + / - 2n] /x dV, 



v 



where dV(£) = \/|g(£)l^£ 1S the volume form, integration is taken over com- 
pact V and i-R is the scalar curvature computed for V. For a real evolution 
parameter r > 0, it is considered J v adV = 1 when fi = (47rr) _ra e~f . 
The functional approach can be redefined for N-anholonomic manifolds 
with stochastic generating functions (in this section, we shall follow the 
N-adapted geometric formalism elaborated in Refs. @[B]): 

Claim 4.1 Definition: For nonholonomic manifolds of even dimension 
In with stochastically generated geometric objects, the stochastic Perelman's 
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Junctionals for the canonical d-connection D are defined 



Hz J) 

W(g,/,r) 



R + S + 



D/ 



~ f dV, 



v 



r [R + S + 



l Df 



+ 



>Df 



+ f-2n 



(26) 
fi dV, (27) 



v 



where dV is the volume form of L g, R and S are respectively the h- and 
v-components of the curvature scalar of D, for D a = (Di,D a ), or D = 



( h D, V D), 



l Df 



+ 



>Df 



an 



d f satisfies J v fidV = 1 f 



or 



fi = (Attt) n e * and r > 0. 



Proof. The formulas (|25p are redefined for some / and / (which can be 
a non-explicit relation between random functions etc and with gravitational 
" noise" induced by stochastic N-adapted R + S) when 



,i?+|V/| 2 e 



-/ 



R + S + 



h Df 



+ 



>Df 



-f 



Re-scaling the evolution parameter r (it is similar to the considered above 
X), t — > r, we have 



r( ,i?+|V/|) 2 + /-2n) 













\ 2 - 


I" = 


T (R + S + 


h Df 


+ 


v Df 


J +/-2n 



for some q and q\ for which J qdV = and j q\dV = 0. □ 

v v 

The geometric objects defining functionals J- and W are with some com- 
ponents computed as expectation values (see details in pQ, for instance, for- 
mulas (12) - (14) when the probability density function subjected to Focker- 
Plank conditions, is used for computing such values). We consider the h- 
variation h 5gij = Vij, the v-variation v 5g a i, = v a b, and h 5f = h f, v 6f = v f 

An explicit calculus for the first N-adapted variations of (j26j) , see similar 
details in [61 [8], distinguished into sure and random components, can be 
represented in the form 



5F( Vij ,v ab , h f, v f) 



{[-v ij {R ij + D i D j f) + {- 



+ [-Vab{Rab + D a D b f) 



(28) 



f)(2 h Af-\ h D f\)+R 



f) (2 v Af-\ v Df\)+S}}e-?dV ; 
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where h A = D l D i and V A = D a D a , A = h A + V A, and h u = g ij v ih v v = 
h ab v ab . 

4.2 Main Theorems for stochastic Ricci flow equations 

We shall prove that stochastic equations of type (|13j) can be derived from 
the Perelman's N-adapted functionals (|26p and (|27p (for simplicity, we shall 
not consider the normalized term and put A = 0). 

Definition 4.1 A metric g generated by a stochastic generating function 
<p(x) is called a (nonholonomic) stochastic breather if for some xi < Xi an d 
a > the metrics a g(xi) and a g(x2) differ only by a N-adapted diffeo- 
morphism. The cases «=,<,> 1 define correspondingly the steady, shrink- 
ing and expanding breathers ( there are configurations when, for instance, the 
h-component of metric is steady but the v-component is shrinking). 

Clearly, the breather properties, in sure and stochastic variables, depend 
on the type of connections which are used for definition of Ricci flows. We 
can elaborate a unique nonlinear diffusion scenarios for D for the assump- 
tions of Hypothesis [TJ 

Following a N-adapted variational calculus for J-"(g, /), with formula 
(|28p . Laplacian A and h- and v-components of the Ricci tensor, Rij and 
Sij, and considering parameter t(x), dr/dx = —1, we formulate 

Theorem 4.1 The stochastic nonholonomic Ricci flows can be parametrized 
by corresponding nonholonomic distributions and characterized by evolution 
equations 



-jij 



dx dx 



op d/tqfe _ g 



dj_ 

dx 



2 



R-S (29) 



-A/ + 

(JX 

and the property that 

q-?( g(x),/(x)) = 2 J [\Ri, + DiDjf\ 2 + \R ab + D a D b f\ 2 \ e~Uv, 



J e 1 dV is constant and the geometric objects T>,A,Rij,R ab ,R and S are 
v 

induced from random generating functions in g(x). 
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Proof. We sketch the idea of such a proof: We should follow G. Perel- 
man [43j constructions (details are given for the connection V in the Propo- 
sition 1.5.3 of [36]) but for the canonical d-connection D). The random 
components of metrics and connections are included as mathematical expec- 
tations. The most important task is to prove that for stochastic generating 
functions the equation (|29p has a unique solution for well defined conditions. 
Really, this equation is equivalent to ({23} for with a stochastic system (f22l) 
can be associated. The existence of a unique solution follows from Theorem 

ED □ 

We can also derive stochastic evolution equations, following stochastic 
N-adapted modifications of Proposition 1.5.8 in |46j containing the details 
of the original result from [13] : 

Theorem 4.2 If a d-metric g(x) determined by stochastic generating func- 
tion (p(x) (Q2P and functions fix) an -d t(x) evolve stochastically following 
the equations ( a stochastic nonholonomic generalization of Hamilton's equa- 
tions) 

d 9ij _ 9 g dhab _ f> 
ox ox 

-L = -A/+ D/ -R-S+-, (30) 
OX r 

I ~ 1 

and the property that that the stochastic Ricci flow evolution is constrained 

to satisfy, for J (47TT)~ n e~^ dV = const, the conditions 
v 

g(x),/(x),r(x)) = 

2 J ?[\R i3 + DiDj- ^ 9ij \ 2 + \R ab + D a D b f- ^g ab \ 2 ](A^)- n e-hv. 
v 

The equation (|30p is similar to (|24p and can be related to a stochastic 
system of type (|22l) . The condition (j3~T1) states a relation between a Ricci 
flow evolution parameter f and stochastic evolution parameter x- 

4.3 Statistical analogy for stochastic Ricci flows 

The Ricci flow theory has, in its turn, a very interesting application in 
the theory of stochastic equations and diffusion. It allows us to characterize 
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additionally such processes via associated statistical functionals, entropy and 
thermodynamical values. G. Perelman emphasized |43j that the functional 
iW is in a sense analogous to minus entropy. In this section, we prove that 
such a property exists also for nonholonomic stochastic Ricci flows and that 
we can provide a statistical model for nonlinear diffusion processesjf] 

We consider a evolution of stochastic geometric systems described by 
some metrics g(r), N-connections N?(t) and related canonical d-connections 
D(t) when the conditions of Theorem 14.21 are satisfied. It follows: 



Theorem 4.3 Any family of stochastically Ricci flow evolving nonholo- 
nomic geometries, and solutions of Einstein equations, with nonlinear dif- 
fusion data derived from Hypothesis [7] and conditions of Theorem 13.11 is 
characterized by thermodynamic values 



(E) = -t z J | // S + 
V 

S = - I ?[ R + S + 



h Df 



+ 



n 



Df + v Df - - U dV, 



v Df 



+ f-2n 



J2 dV, 



(32) 



d = 2r 4 



v 



\Rij + DiDjf- — gij \ 2 + \R ab + D a D b f- ^g ab \ 



Ji dV. 



Proof. It follows from a straightforward computation for 
Z = exp{J v [— / + n]fidV} as in the original paper [43]. For nonholonomic 
stochastic processes, we have to N-adapt the constructions as in [U[8]. The 
stochastic terms are included in formulas (|32|) via expected values of smooth 
coefficients satisfying the Fokker-Plank equation (or forward Kolmogorov 
equation). □ 

Any N-adapted stochastic configuration determined by a canonical d— 
connection D is thermodynamically more (less, equivalent) convenient than 
a similar one defined by the Levi-Civita connection V if S < ,S (S > 
,S,S = ,S). Similarly, a certain geometry can be more (less, equivalent) 
convenient than a stochastic one and related Ricci flow evolution models. 



Let us remember some important concepts from statistical mechanics: The partition 
function Z — f exp(—/3E)du)(E) for the canonical ensemble at temperature T /3 _1 is 
defined by the measure taken to be the density of states uj(E). The thermodynamical 
values are computed in the form: the average energy, (E) — —d\ogZ/d T f3, the entropy 
S = T /3(E) + log Z and the fluctuation a = {(E - (E)) 2 ) = d 2 log Z/d T /3 2 . 
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